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ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ 
В МОДЕЛИ ЭНДОГЕННОГО РОСТА 
ЗАМКНУТЫХ ЭКОНОМИЧЕСКИХ СИСТЕМ 


В статье рассматривается математическая модель управления распределением ин- 
вестиций, в которой научно-технический прогресс учитывается в форме увеличе- 
ния коэффициента эластичности по фактору капитал. Предложен алгоритм реше- 
ния соответствующей задачи управления, при котором средства выделяются из 
объёма, выделявшегося ранее на инвестирование в производство. В результате в 
конечный момент времени удаётся достичь увеличения объёма выпуска по сравне- 
нию с объёмом, полученным по итогам рассматриваемого периода, а также полу- 
чить максимальное накопление капитала. 

Ключевые слова: оптимальное управление, критерий максимальности. 


Введение. В работах Р.Солоу показано, что для неоклассических произ- 
водственных функций экономический рост в долгосрочном плане ограни- 
чен темпами прироста трудовых ресурсов [1]. Поэтому, вскоре, внимание 
экономистов, занимающихся данной проблематикой, было переключено на 
научно-технический прогресс (НТП) как основной источник экономического 
роста. Представители нового направления неоклассицизма сводили НТП к 
экзогенному феномену. В 50-е годы 20-го столетия появились модели учёта 
НТП в виде введения в производственную функцию некоторого параметра 
А(Е). В зависимости от места привязки коэффициента А стали различать 
«нейтральный», по Хиксу, Харроду и Солоу, НТП [2]. Но при такой трактов- 
ке вопроса оставались в стороне проблемы мотивации инновационного 
процесса, формирования параметра А, и сам этот параметр являлся возрас- 
тающей функцией времени. Желание разрешить поставленные проблемы 
привело к отказу от концепции экзогенного НТП. 

Теория эндогенного экономического роста появилась в конце 
80-х гг. 20-го века. Переходными моделями можно назвать модели Эрроу 
[3] и Шешински [4], в которых НТП учитывался в процессе обучения при 
использовании основных фондов или же с помощью некоторого фактора А, 
определяемого уровнем знаний. В моделях Лукаса [5], Ромера [6,7], Агито- 
на-Ховитта [8] влияние инноваций на технологические сдвиги моделиро- 
валось с учётом накопления человеческого капитала с помощью блока про- 
межуточных товаров, влияющего на качество конечного продукта, себесто- 
имость производства, эффективность основного капитала. Этот же подход 
развивается в работах Узавы [9], Беккера [10] и других представителей 
этого направления. Более подробный обзор этих моделей содержится, 
например, в [11]. 

Несмотря на разнообразие и многочисленность предложенных под- 
ходов, исследование внутренних механизмов влияния научно-технического 
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прогресса на эндогенный экономический рост всё ещё остаётся одной из 
главных нерешённых проблем в моделировании экономических систем. 
Другими словами, до сих пор не раскрыта суть той движущей силы, той 
«пружины», которая, подобно пружине часового механизма, является пер- 
вопричиной постоянного экономического роста и социального прогресса. В 
работе авторов [12] показывается, что рост эластичности производствен- 
ной функции по фактору капитал в рамках неоклассической модели роста 
Солоу может приводить к существенному росту накоплений. Отталкиваясь 
от предложенного в [12] подхода, в данной работе предлагается модель 
оптимального управления, в которой можно добиться увеличения выпуска 
и накопления капитала, выбирая некоторую последовательность управле- 
ний, направленных на получение эффективных технологий. Формулируется 
некоторая задача оптимизации, решение которой позволяет установить оп- 
тимальный уровень инвестиций, необходимый для достижения максималь- 
ного темпа экономического роста. Данная модель тестируется на примере 
советской экономики за период 1961-1985 гг. Показано, что увеличение ин- 
вестиций в научно-исследовательский сектор (НИС) могло бы привести к 
увеличению темпов роста на 10-15 %. 

Постановка задачи. Из результатов, изложенных в [12], следует, 
что при Ч х > 1, {> ® функция накопления капитала К(Е), удовлетворя- 
ющая уравнению Рамсея — Солоу, теоретически может расти неограничен- 
но (при этом соответствующая производственная функция удовлетворяет 
всем условиям неоклассической функции). Легко показать, что вообще с 


любым ростом накопление К(Ё) возрастает. Поэтому одна из возмож- 
ных трактовок эндогенного экономического роста может быть связана с 
тем, что научно-технический прогресс приводит к росту | '‹ |, т.е. росту 
эластичности по фактору капитал, другими словами, приводит к повыше- 
нию эффективности отдачи в использовании основных фондов. Тогда ясно, 
что инвестиции в исследования и новые разработки в научно-исследова- 
тельский сектор будут приводить к росту | ‹ | и, как следствие, к росту 
К(Е. Это и есть один из важнейших механизмов эндогенного экономическо- 
го роста. В такой постановке возникает естественная задача оптимального 
управления параметром | ‹ 

Будем считать, что для абстрактного производственного цикла 
можно получить некоторую функцию, аппроксимирующую значения коэф- 
фициента эластичности по фактору капитал | ** | от расходов на науку °„, 
т.е. от объёма инвестиций $, Е где Е - объём выпуска в предшествующий 
момент времени. Тогда можно считать, что у нас есть некоторая 
возрастающая функция @ к =@ к (5. Е). Таким образом можно получить 
для каждого управления $„ некоторую качественную оценку изменения ко- 
эффициента эластичности по фактору капитаа ЯАк и рассмотреть 
производственный цикл того же продукта, но уже с другими технологиче- 
скими условиями. 
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Заметим, что в рассматриваемой модели оптимального управления 
фактор труда Г.(#) считается некоторой заранее известной функцией, изме- 
нение которой в большой замкнутой системе определяется в основном де- 
мографическим фактором. 

Пусть при :<к имеем некоторый первоначальный производственный 
процесс. В результате накопления на этом этапе фактор капитал достиг 
значения Ко. Пусть при этом объём трудовых ресурсов равен Г. Тогда на 
следующем временном интервале /1,;1/] система развивается согласно 
уравнению Солоу [1], в котором для определенности принимается произ- 
водственная функция Кобба-Дугласа. Пусть при этом 5 — доля выпуска, ин- 
вестируемая в старые основные фонды (будем считать ее константой). С 
целью повышения эффективности и достижения более оптимального 
управления будем считать, что величина 5,= 5, 0< 5 <5 обозначает инве- 
стиции в научно-исследовательский сектор и изымается из средств, ранее 
инвестировавшихся в старые основные фонды. Тогда согласно модели Со- 
лоу имеем нелинейное обыкновенное дифференциальное уравнение, опи- 
сывающее динамику капитала с течением времени, 

= (5-5) АКП -ц К(Ю, Г= И, (1) 
и 
которое мы рассматриваем на временном отрезке &<1<#. Здесь А>0 - не- 
который известный числовой коэффициент; д - коэффициент амортизации; 
0<0<1 — коэффициент эластичности по фактору капитал. 

Если данный отрезок времени мал по сравнению со всем периодом 
цикла производственного развития, то все функции, входящие в правую 
часть уравнения Солоу (1) на данном временном интервале, можно считать 
постоянными. Можно доказать, что при этом погрешность в определении 
приращения капитала является малой второго порядка. При таких условиях 
можно считать функцию [() имеющей некоторое постоянное значение [.. 

Уравнение (1) с начальным условием 


К(к)=К, (2) 
представляет собой классическую задачу Коши, решение которой имеет 
следующий вид: 


К(Ё)= (5- а ме 


1- е = №-0] + К (1-в и - 1) „Ч =во. (3) 
д 


Из этого явного вида легко извлекается новое значение капитала 
К, в момент времени #„. Соответственно, легко вычисляется новое значение 
выпуска /Р`в момент времени #,, т.е. 


а (4) 

Зададим новое значение 5,= 5, 0< 5, <1 для промежутка /1,;5]. Из 

функции бк =@Як (5,ЁЕ:) находим новое значение @1. Тогда произ- 

водственный процесс для промежутка #,<#<# может быть описан аналогич- 
но уравнению (1): 
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ак [04 -@\ 
—-= 5-5) АК "И -и КО (5) 
с новым начальным условием 
ЖЕ К . (6) 


Решая новую задачу Коши (5) и (6), получаем новое значение капи- 
тала К› в момент времени 2›, что приводит к новому значению выпуска: 


Е, = ЯК. 1 в в ‚ где Г - объём трудовых ресурсов в момент времени 


р. 

Продолжая процесс и задавая новое значение $» получаем новое 
значение а. и можем решать уравнение, аналогичное (5), для <<. Та- 
ким образом, получим новые значения К; и Р; - соответственно объёмы 
основных фондов и выпуска в момент времени &. 

В общем виде приходим к следующему алгоритму вычисления зна- 
чений капитала и объёма выпуска в последующий момент времени #1»; в 
зависимости от известных значений тех же величин в момент времени {. 
Пусть в момент времени & известны величины: 


К(&)= К, ВЫ НЕЕЕ. (7) 
Фиксируем значение величины 5,=5, т.е. долю выпуска, инвестиру- 
емую в научно-исследовательский сектор, причём 0<5;<5. Таким образом, 
получаем величину 5:2; — некоторый объём инвестиций на науку, который 
позволяет увеличить эффективность капитала, т.е. повысить коэффициент 
эластичности по фактору капитал до значения @ ; =@ к (5, Е). Тогда для 
&<1<., получим следующее дифференциальное уравнение: 
К АК КО. (8) 
И 

Решая задачу Коши, состоящую из уравнения (8) и начального 

условия (7), получаем К’'.,=К(&.-,) и соответственно можем вычислить 
Е = АК!’ р" › 
где [,,., — объём трудовых ресурсов при #=#.... 

Очевидно, что при фиксированных значениях правых частей (7) ре- 
шение (8) зависит только от величины 5. Таким образом, можно, меняя 
значение 5; и решая задачу (7), (8), получать различные значения капита- 
ла (или выпуска) в момент времени #,... 

Отметим, что величина значений К, Ё; напрямую зависит от выбора 
величины 5; на соответствующем временном промежутке /1, 1,.//. Естествен- 
но, что значения итоговых величин К, РЁ, зависят от всего набора значе- 
НИЙ 5, 1=0,...,п-1. Тогда можно поставить следующую задачу: как выбрать 
оптимальное управление распределением инвестиций в НИС, т.е. последо- 


| я ы 
вательность а так, чтобы получить максимальныи экономическии 


эффект, т.е. максимальную величину К). 
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Сформулируем оптимизационную задачу более конкретно: 
К(1)= К, Е )=Ь,Е(1)=Е 


1 


и; =ак (5, Е) 
С - (5- 5) АК -ь К(1), 1 = 0, ...,п-1 (9) 
тах К, - ? 


Здесь КР, Г.,...„.„ — статистические данные, К; для 7>0 — решение (1-1)-й 
задачи Коши, Р; вычисляется по формуле Е = АК; "! и для 1>0. 
Методы решения оптимизационной задачи. Полученную оптимизаци- 


онную задачу можно решать численно различными способами, например, 
методом покоординатного спуска или методом градиентного спуска. В ре- 


зультате решения получаем некоторую последовательность [о т.е. 


получаем план оптимального управления научными исследованиями и раз- 
работками на период времени /1, &,/. При этом критерием выбора элемен- 
тов этой последовательности является максимальность величины основных 
фондов в момент времени #,. Нужно отметить, что в качестве & можно 
рассмотреть годы, кварталы, месяцы в зависимости от объекта исследова- 
ния и наличия необходимых статистических данных. 

Протестируем построенную модель на примере статистических дан- 
ных для СССР за 1961-1985 гг. Пусть известны: К, 1=0,...,п-1 — основные 
фонды за каждый год (в млрд. руб.); Г, 1=0,...,п-1 — объем трудовых ресур- 
сов (пересчитанный через среднюю зарплату в млрд. руб.); 
Е, 1=0,....п-1 — валовой общественный продукт (в млрд. руб.); Г, 1=0,...,и-1 
— расходы на науку из госбюджета и других источников (в млрд. руб.), 
п=25. Эти данные приведены в табл.1. 

Таблица 1 
Основные макроэкономические показатели советской экономики 


Валовой Основные Трудовые Коэффициент Расходы 
продукт Е фонды К ресурсы Ё эластичности на науку \ 
(в млрд.руб.) | (в млрд. руб.) | (в млрд. руб.)| по фактору (в млрд. руб.) 
капитал а 
5 


0,342857 

0,36 

0,375 
0,388235 
04 9 


ен 
0,323077 
[66 | 

7,6 





Окончание табл.1 
5 Г. 20 1 
1968 754,6 114,08 0,410526 





1969 _| 622.2 813,5 123,72 0,42 


114 


Раздел «Управление, вычислительная техника и информатика» 








1970 _| 643,5 860 132,05 0.428571 117 





Считаем, что коэффициент эластичности о меняется за указанный 
период с 0,3 до 0,5. Поэтому можно составить функцию, аппроксимирую- 
щую зависимость коэффициента эластичности по фактору капитал от 
объёма инвестиций в научно-исследовательский сектор, 

_ 0,6+7, 0,615 
ь 7+ 6,95 
Полученную функцию (10) будем использовать для оценки значений «х при 
новых (изменённых) значениях /;=5,Ё; 

Из анализа статистических данных за указанный временной проме- 
жуток получаем, что доля выпуска, инвестируемого в старые основные 
фонды, равна 5=0,165. При этом величина инвестиций в научно-исследова- 
тельский сектор не превышала 1,5% от валового продукта, т.е. х,=0,015. 
Поставим следующую задачу управления: какую часть 5 нужно было инве- 
стировать ежегодно в научно-исследовательский сектор для того, чтобы 
добиться наилучшего экономического эффекта? 

Имеем: К,=441 млрд. руб.; Р.=324 млрд. руб.; [,=65,97 млрд. руб. 
Тогда при 1=0 можно выбрать величину 5%, 5% :0,015< 5, < 0,165. Далее вы- 


(10) 


числяем И = 55 Ро, ао находим по формуле (10) при :=0. В итоге получаем 
следующую задачу Коши: 
ее. (5-5) А Ко 0 -у К 
ий (11) 
К(= 0) = Ку, [(!= 0)= Ц, Е(1= 0)= Е, 
где и=0,05. 

Решая дифференциальное уравнение (11) описанным выше спосо- 
бом, получаем некоторое значение К; . Для новых значений К’ и си можно 
вычислить новое значение Ё = А К° р ‚ Т.е. мы получили новые зна- 
чения основных фондов и валового продукта для второго года из рассмат- 
риваемого временного промежутка, иначе говоря, при 1=1. Теперь обозна- 
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чим Ро=А, Ко=Кь шл=Ё, и, выбирая следующее 5, 1=1, считаем а, по (10) 
и решаем новую задачу (5), (6). Проводим аналогичные рассуждения для 
1=2,...п-1. В нашем случае и=25. Таким образом, мы можем составить по- 


5 я в 
следовательность |5; — и для каждои такои последовательности получа- 
1/1=0 


ем свои значения К,, Р‚. 

Отметим, что на основе имеющихся статистических данных для ука- 
занного временного промежутка удалось подобрать функцию, аппроксими- 
рующую величину Г в зависимости от времени. Эта функция не является 
константой и имеет следующий вид: 


[(1= 61+ 51+ 0,18Г - 0,0015Ё. (12) 


Используем соотношение (12) для подсчёта Г,, 1=0,...,п-/. 

Отметим, что для лучшей аппроксимации статистических данных 
величину 4 также следует считать переменной. Нам удалось подобрать для 
нее следующее аппроксимирующее выражение: 

А= 2,99- 0,0363, 


где 1=0,...,и-1. Легко показать, что рассуждения, сделанные выше, будут 
верны и в этом случае. 

При сделанных выше пояснениях задача (9) оптимального управле- 
ния была решена двумя разными методами: методом покоординатного 
спуска (МКС) и методом градиентного спуска (МГС). Из двух указанных ме- 
тодов более эффективным оказался МГС. Следует заметить, что, поскольку 
целевой функционал не удается выразить в явном виде, для реализации 
МГС значения градиента в текущей точке находили численно по формуле 
конечных разностей. 

Для дополнительного тестирования эффективности предложенных 
алгоритмов был применен один из вариантов глобального случайного по- 
иска [13] при решении сформулированной выше оптимизационной задачи. 
Основанием для этого послужило то, что описанные выше два метода но- 
сят локальный характер и в оптимизационных задачах со сложным поведе- 
нием целевой функции могут приводить к нахождению лишь локальных 
максимумов, которые часто не совпадают с глобальным. Алгоритм, описан- 
ный в [13], представляет собой пошаговый итерационный процесс. На каж- 
дом шаге при нахождении максимума многомерной целевой функции Р\(х), 
которая для простоты считается положительной, происходит вычисление 
ее значений в М случайных точках: ЕР, = Е(х„),т = 1....,М. Далее подсчи- 


тываются вероятности 





бы 
Ри = ‚т= 1... М. 


На следующем шаге итерации число новых случайных точек в 
окрестности «хороших» пробных точек предыдущих итераций выбирается 
пропорционально этим вероятностям. Таким образом, пробные точки сле- 
дующих шагов выбираются в окрестности самых лучших пробных точек 
предыдущих шагов. Кроме того, области, в которых выбираются случайные 
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точки последующих итераций, постепенно стягиваются к малым окрестно- 
стям «хороших» точек предыдущих итераций. 

Реализация данного алгоритма на рассматриваемой тестовой зада- 
че из экономики СССР привела к тем же числовым результатам, что и два 
вышеописанных локальных метода. 


Результаты численного моделирования. В табл.2 указаны реальные 
п-1 
1=0 
методом градиентного спуска. 


значения (5,) в 1961-1985 гг и значения тех же величин, полученных 


Таблица 2 
Зависимость основного капитала от инвестиционной политики 


5, Основные фонды, |9, 1=0,,...,П-1, Основные фонды, 
1=0,п-1, в млрд. руб., полученные мето- | в млрд. руб., полу- 
по стат. по стат.данным дом градиентного | ченные методом 
данным спуска градиентного 


спуска 
0,015 442 
0,015 477 
0,015 515 
0,015 557 
0,015 601 
0,015 649 
0,015 700 
0,015 757 
0,015 814 
0,015 860 [Г 0020 [|965 | 
0,015 914 
0,015 972 
0,015 1082 
0,015 1150 
0,015 1256 
0,015 1345 
0,015 1437 
0,015 1537 
0,015 1638 
0,015 1742 
0,015 1851 
0,015 1968 
0,015 2092 
0,015 2234 
1985 [0,015 2333 





Из табл.2 видно, что решение поставленной оптимизационной зада- 
чи позволяет получить значительное увеличение целевой функции. Реали- 
зация данного алгоритма на рассматриваемой тестовой задаче из экономи- 
п-1 
1=0' 
рых имеем К›;=2576 млрд. руб. (этот результат получен во всех трех мето- 
дах). Таким образом, управляя инвестиционным процессом и выделяя на 
научные исследования и разработки указанные доли валового продукта, 


ки СССР привела к получению параметров управления (5;) при кото- 
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можно было бы получить значение на 10,4% превышающее реальное зна- 
чение капитала в 1985 году, которое реально составило 2333 млрд. руб. 
Выводы. Предложенная математическая модель основана на предположе- 
нии о том, что более эффективные технологии приводят к увеличению та- 
кого макроэкономического параметра, как коэффициент эластичности по 
фактору капитал. При этом можно подобрать такую последовательность 
управлений, которая максимизирует накопление капитала. Предложенный 
алгоритм может быть использован для решения задач оптимизации в си- 
стемах, описывающихся с помощью дифференциальных уравнений, опре- 
деляющих поведение системы во времени. При этом от решения диффе- 
ренциального уравнения на предыдущем временном интервале зависят на- 
чальные условия для задачи Коши на последующем интервале. Если это 
решение зависит от некоторых параметров управления, то предложенный 
алгоритм позволяет найти такую последовательность этих параметров, при 
которой достигает максимума некоторый критерий, выбранный в качестве 
оценки эффективности данной системы. 

Автор выражает благодарность профессору М.А. Сумбатяну за по- 
становку задачи. 
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ОРТТМАЕ СОМТКОЕ ТМ А МОБЕЁ ОЕ ЕМООСЕМОЦ$ СКОМ/ТН 
РОК СГО$ЕО ЕСОМОМТС 5У5ТЕМ$ 
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та еа$ЧсКу соейаепе оп Ше саркКа! Гас{ог. ТПеге 15 ргорозе а зо]иЧоп а!до- 
Гбт {0 эомМе гезресИиуе орЧтмаНоп сопёго! ргоБет, м/Пеге {Ле шуезётепе 
Рипаб аге {аКеп гот {Не аточцпЕ м/МСП ма еа ег и{епае4 Гог пуезтет © 
{Не тапуРасиге ргосез$. Аз а гезий, опе сап асШеуе а{ а Ипа! тотепё оЁ 
Чте ап ад4 юпа! ходите оЁ ргодисйоп, мПеп сотрагед мА {1е опе д№еп Бу 
геа| 5аНзЫса| Ча{а, аз ме! аз {© геасй а тахитит сарКа! асситцайоп. 
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